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Abstract 

We report the proof that the expression of extended Gibrat's law is unique and the 
probability distribution function (pdf) is also uniquely derived from the law of detailed 
balance and the extended Gibrat's law. In the proof, two approximations are employed that 
the pdf of growth rate is described as tent-shaped exponential functions and that the value 
of the origin of growth rate is constant. These approximations are confirmed in profits data 
of Japanese companies 2003 and 2004. The resultant profits pdf fits with the empirical data 
with high accuracy. This guarantees the validity of the approximations. 
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1 Introduction 

In the large scale region of income, profits, assets, sales and etc (x), the cumulative probability 
distribution function (pdf) P{> x) obeys a power-law for x which is larger than a certain 
threshold Xq: 

P(> x) oc x^V" for X > Xq . (1) 

This power-law and the exponent /i are called Pareto's law and Pareto index, respectively [1]. 
The power- law distribution is well investigated by using various models in econophysics [2] . 

Recently, Fujiwara et al. [3] find that Pareto's law can be derived kinematically from the law 
of detailed balance and Gibrat's law [4] which are also observed in the large scale region x > xq. 
In the proof, they assume no model and only use these two laws in empirical data. 

The detailed balance is time-reversal symmetry (x\ «-> x?): 

P\2{xi,x 2 ) = Pvz(x 2 ,xi) . (2) 



* e-mail address: ishikawa@kanazawa-gu.ac.jp 



1 



Here x\ and X2 are two successive incomes, profits, assets, sales, etc. and P±2(xi, x 2 ) is the joint 
pdf. Gibrat's law states that the conditional pdf of growth rate Q(R\x\) is independent of the 
initial value X\\ 

Q(R\ Xl ) = Q(R) . (3) 

Here growth rate R is defined as the ratio R = X2/X1 and Q(R\x±) is defined by using the pdf 
P(xi) and the joint pdf Pir(xi, R) as Q(R\xi) = Pir{x\,R)/P(xi). 

In Ref. [5], the kinematics is extended to dynamics by analyzing data on the assessed value 
of land in Japan. In the non-equilibrium system we propose an extension of the detailed balance 
(detailed quasi-balance) as follows 

Pi2(xi,x 2 )=p 12 ((^y 9 , axi ) . (4) 

From Gibrat's law (3) and the detailed quasi-balance (4), we derive Pareto's law with annually 
varying Pareto index. The parameters 6, a are related to the change of Pareto index and the 
relation is confirmed in the empirical data nicely. 

These findings are important for the progress of econophysics. Above derivations are, how- 
ever, valid only in the large scale region where Gibrat's law (3) holds. It is well known that 
Pareto's law is not observed below the threshold xo [4, 6]. The reason is thought to be the 
breakdown of Gibrat's law [3]. The breakdown of Gibrat's law in empirical data is reported by 
Stanley's group [7]. Takayasu et al. [8] and Aoyama et al. [9] also report that Gibrat's law does 
not hold in the middle scale region by using data of Japanese companies. 

In Ref. [10], Gibrat's law is extended in the middle scale region by employing profits data 
of Japanese companies in 2002 and 2003. We approximate the conditional pdf of profits growth 
rate as so-called tent-shaped exponential functions 

Q(R\ Xl ) = d(rci) iT^O-i for R > 1 , (5) 
Q(R\ Xl ) = cf(xi) iZ+M*!)- 1 f or R < 1 . ( 6 ) 

By measuring t± we have assumed the x\ dependence to be 

X\ 

t±(xi) = t±(x ) ±a± In — , (7) 

Xo 

and have estimated the parameters as [11] 

a + ~ a_ ~ for x\ > xo , (8) 

q + ~ a_ 7^ for x m i n < xi < xo , (9) 

t+(xo)-Mx ) ~ M- (10) 

From the detailed balance (2) and extended Gibrat's law (7) - (10), we have derived the pdf in 
the large and middle scale region uniformly as follows 

P{x) = Cx-^ 1 ) e ' aln2 1 for x > x min , (11) 
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where a = (a + + a_) /2. This is confirmed in the empirical data. 

In this study, we prove that the x\ dependence of t± (7) with a + = q_ is unique if the 
pdf of growth rate is approximated by tent-shaped exponential functions (5), (6). This means, 
consequently, that the pdf in the large and middle scale region (11) is also unique if the x\ 
dependence of d(x\) is negligible. We confirm these approximations in profits data of Japanese 
companies 2003 and 2004 [12] and show that the pdf (11) fits with empirical data nicely by the 
refined data analysis. 

2 Growth rate distributions of profits in the database 

In the database, Pareto's law (1) is observed in the large scale region whereas it fails in the 
middle one (Fig. 1). At the same time, it is confirmed that the detailed balance (2) holds not 
only in the large scale region x\ > xq and xi > xq but also in all regions x\ > and xi > 
(Fig. 2). 1 

The breakdown of Pareto's law is thought to be caused by the breakdown of Gibrat's 
law in the middle scale region. We examine, therefore, the pdf of profits growth rate in 
the database. In the analysis, we divide the range of x\ into logarithmically equal bins as 
xi € 4 x [loi+o-^™- 1 ), l0 1+0 - 2ri ] thousand yen with n = 1, 2, ■ ■ ■ , 20. In Fig. 3-6, the probability 
densities for r are expressed in the case of n = 1, • • • , 5, n = 6, • • • , 10, n = 11, • • • , 15 and 
n = 16, • • • , 20, respectively. The number of the companies in Fig. 3 - 6 is "22, 005", "89, 507", 
"85,020" and "24,203", respectively. Here we use the log profits growth rate r = log 10 R- The 
probability density for r defined by q{r\x\) is related to that for R by 

log 10 Q(R\ Xl ) + r + log 10 (ln 10) = log 10 q(r\ Xl ) . (12) 

From Fig. 3-6, log 10 q{r\x\) is approximated by linear functions of r as follows 

l°gio l{ r \ x i) = c ( x i) — t+( x i) r for r>0, (13) 
log 10 g(r|xi) = c(xi) + t-(xi) r for r<0. (14) 

These are expressed as tent-shaped exponential functions (5), (6) by d{x\) = lO^^/ln 10 . In 
addition, the x\ dependence of c{x\) (d(xi)) is negligible for n = 6, • • • , 20. The validity of these 
approximations should be checked against the results. 

3 The uniqueness of extended Gibrat's law 

In this section, we show that the x\ dependence of t± (7) is unique under approximations 
(5), (6) ((13), (14)). 

Due to the relation of ^12(^1, X2)dx\dx2 = P\r{x\, R)dx\dR under the change of variables 
from (xi, X2) to (xi,R), these two joint pdfs are related to each other Pm(xi,R) = xiPu(xi,X2)- 

lr The scatter plot in Fig. 2 is different from one in Ref [10]. The reason is that the identification of profits in 
2002 and 2003 in Ref. [10] was partly failed. As a result, the pdfs of profits growth rate are slightly different from 
those in this paper. The conclusion in Ref. [10] is, however, not changed. 
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By the use of this relation, the detailed balance (2) is rewritten in terms of Pir(x±, R) as follows: 

P 1R { Xl ,R) = R- 1 P 1R {x 2 ,R- 1 ). (15) 

Substituting the joint pdf P\r(xi,R) for the conditional probability Q{R\x\), the detailed bal- 
ance is expressed as 

P( Xl ) 1 Q{R- l \x 2 ) 



P(x 2 ) R Q(R\xi) ' 
Under approximations (5) and (6), the detailed balance is reduced to 



P(xi) d(x 2 ) 



U+t+(xi)-t_(x2)+l 



P{x 2 ) d{ Xl ) 

for R > 1. By using the notation P(x) = P(x)d(x), the detailed balance becomes 

p ( x ±) = R +t+{ Xl )-t-{R xi)+l 

P(R Xl ) 



(16) 



(17) 



(18) 



By expanding Eq. (18) around R = 1, the following differential equation is obtained 

l + t+(x) - t_(x)]P(x) +x P\x) = 0, (19) 

where x denotes x\. The same differential equation is obtained for R < 1. The solution is given 
by 



P(x) = Cx- 1 e~ G{x) , 



(20) 



where t + (x) — t-(x) = g(x) and J g(x)/x dx = G(x). 

In order to make the solution (20) around R = 1 satisfies Eq. (18), the following equation 
must be valid for all R: 



-G{x) + G{Rx)= t+(x)-t-(Rx) 
The derivative of Eq. (21) with respect to x is 



lni? . 



g{x) g{R x) 



t+ (x) - Rt- (R x) 



InR . 



x x 

By expanding Eq. (22) around R = 1, following differential equations are obtained 

x t+"(x) +t-"(x) +t+'(x) +t-'(x) = , 
2 t+'(x) + tJ(x) - 3x tJ\x) - x 2 \t + (3) (x) + 2 tJ 3) (x) 

The solutions are given by 

t+(x) = -^ln 2 z + (C + i-C_i)lnx+(C + o-C_o) 
C 

t_( x ) = — ^ln 2 x + C_ilna; + C_o • 



(21) 

(22) 

(23) 
(24) 

(25) 
(26) 



To make these solutions satisfy Eq. (21), the coefficients must be C- 2 = and C+i = 0. Finally 
we conclude that t±(x) is uniquely expressed as Eq. (7) with a + = a_. 



4 The profits distribution and the data fitting 



Under approximations (5) and (6) ((13) and (14)), we obtain the profits pdf 

P(x) = P{x)d{x) = Cx-^ + V e ~ aln2 =S , (27) 

where we use the relation (10) confirmed in Ref. [11]. 

In Fig. 7, t± hardly responds to x\ for n = 17, • • • , 20. This means that Gibrat's law holds 
in the large profits region. On the other hand, t + linearly increases and t_ linearly decreases 
symmetrically with log 10 x\ for n = 9, 10, • • ■ , 13. The parameters are estimated as Eq. (8) and 
(9) with a (= a + = a_) ~ 0.14 and x = 4 x 10 1+a2 ( 17_1 ) ~ 63, 000 thousand yen. Because the 
x\ dependence of c{x\) (d{x\)) is negligible in this region, the profits pdf is reduced to Eq. (11). 
We observe that this pdf fits with the empirical data nicely in Fig. 8. 

Notice that the estimation of a in Fig. 7 is significant. If we take a slightly different a, the 
pdf (11) cannot fit with the empirical data (a = 0.10 or a = 0.20 in Fig. 8 for instance). 

5 Conclusion 

In this paper, we have shown the proof that the expression of extended Gibrat's law is unique 
and the pdf in the large and middle scale region is also uniquely derived from the law of detailed 
balance and the extended Gibrat's law. In the proof, we have employed two approximations that 
the pdf of growth rate is described as tent-shaped exponential functions and that the value of 
the origin of growth rate is constant. These approximations have been confirmed in profits data 
of Japanese companies 2003 and 2004. The resultant pdf of profits has fitted with the empirical 
data with high accuracy. This guarantees the validity of the approximations. 

For profits data we have used, the distribution is power in the large scale region and log- 
normal type in the middle one. This does not claim that all the distributions in the middle 
scale region are log-normal types. For instance, the pdf of personal income growth rate or sales 
of company is different from tent-shaped exponential functions [3]. In this case, the extended 
Gibrat's law takes a different form. In addition, we describe no pdf in the small scale region [13]. 
Because the x\ dependence of d(x\) in this region is not negligible (Fig. 3). 

Against these restrictions, the proof and the method in this paper is significant for the 
investigation of distributions in the middle and small scale region. We will report the study 
about these issues in the near future. 
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Figure 1: Cumulative probability distributions P(> x\) and P(> X2) for companies, the profits 
of which in 2003 (x\) and 2004 (X2) exceeded 0, x\ > and X2 > 0. The number of the 
companies is "227,132". 
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Figure 2: The scatter plot of all companies in the database, the profits of which in 2003 (xi) 
and 2004 (^2) exceeded 0, x\ > and X2 > 0. The number of the companies is "227,132". 




Figure 3: The probability density distribution q{r\x\) of the log profits growth rate r = 
log 10 x-ijx\ from 2003 to 2004. The data points are classified into five bins of the initial profits 
with equal magnitude in logarithmic scale, x\ € 4 x [io 1 +°- 2 ( n ^ 1 ) ) 10 1+0 - 2 ™] (n = 1,2, •••,5) 
thousand yen. The number of companies in this regime is "22,005". 
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Figure 4: The probability density distribution q{r\x{) of the log profits growth rate r = 
log 10 X2/X1 from 2003 to 2004. The data points are also classified into five bins of the initial prof- 
its with equal magnitude in logarithmic scale, x\ £ 4 x [I0 1+a2(n_1 ), io 2+0 - 2ri ] (n = 6, 7, ■ ■ ■ , 10) 
thousand yen. The number of companies in this regime is "89,507" . 
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Figure 5: The probability density distribution q{r\x\) of the log profits growth rate r = 
log 10 X2/X1 from 2003 to 2004. The data points are also classified into five bins of the initial prof- 
its with equal mag nitude in logarithmic scale, x x € 4x [io^ - 2 ^- 1 ) , io 1+0 - 2 ™] (n = 11, 12, • • • , 15) 
thousand yen. The number of companies in this regime is "85,020". 
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TheLogarithmofProfitsGrowthRater=log 
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Figure 6: The probability density distribution q{r\x\) of the log profits growth rate r = 
log 10 X2/X1 from 2003 to 2004. The data points are also classified into five bins of the initial prof- 
its with equal mag nitude in logarithmic scale, x x € 4x [lO 1 * - 2 ^ -1 ), 10 1+0 - 2 ™] (n = 16, 17, • • • , 20) 
thousand yen. The number of companies in this regime is "24,203" . 
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Figure 7: The relation between the lower bound of each bin x\ and t±(x±). From the left, each 
data point represents n = 1, 2, • • • , 20. 
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Figure 8: The probability distribution function (pdf) P{x2) for companies, the profits of which 
in 2003 (xi) and 2004 (xz) exceeded 0, X\ > and xi > 0. The pdf derived from the detailed 
balance and the extended Gibrat's law fits with the data accurately. Indices /j,, a and xq 
are already estimated in the extended Gibrat's law. If the parameters are different from the 
estimation, the pdf cannot fit with the data (a = 0.10 or a = 0.20 for instance). 
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